
alpha = -180:5:180;
beta = -90:5:90;

alpha_r = deg2rad(alpha);
beta_r = deg2rad(beta);

[alpha_mesh, beta_mesh] = meshgrid(alpha_r, beta_r);

R1 = abs(cos(alpha_mesh).*1);

alphaf=deg2rad(90);
betaf=deg2rad(0);

R=[1/4 -1/4 0; %матрица координат элементов
    -1/4 -1/4 0;
    -1/4 1/4 0;
     1/4 1/4 0;];
k=[0 cos(alphaf) sin(alphaf)];%источник

dphi=2*pi*k*R';%набег фаз

B=exp(-1i*dphi);%вектор фокусировки

for i=1:size(alpha_mesh,1)
for j=1:size(beta_mesh,2)
alpha=alpha_mesh(i,j);
beta=beta_mesh(i,j);
k1=[cos(alpha)*sin(beta) sin(alpha)*sin(beta) cos(beta)];
dphi1=2*pi*k1*R';
B1=exp(1i*dphi1);
F(i,j)=(abs(sum(B.*B1)))^2;
end;
end;

[x, y, z] = sph2cart(beta_mesh, alpha_mesh, F);

figure(1);
surf( x, y, z)
xlabel('x'); ylabel('y'); zlabel('z')

minc = min( [min(min(x)) min(min(y)) min(min(z))] );
maxc = max( [max(max(x)) max(max(y)) max(max(z))] );
xlim([minc maxc]); ylim([minc maxc]); zlim([minc maxc]);

% surf( alpha_mesh, beta_mesh, R)